Implementation of holonomic quantum gates by an isospectral deformation 

of an Ising dimer chain 
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We exactly construct one- and two-qubit holonomic quantum gates in terms of isospectral de- 
formations of an Ising model Hamiltonian. A single logical qubit is constructed out of two spin-^ 
particles; the qubit is a dimer. We find that the holonomic gates obtained are discrete but dense 
in the unitary group. Therefore an approximate gate for a desired one can be constructed with 
arbitrary accuracy. 
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I. INTRODUCTION 

A reliable implementation of a quantum gate is re- 
quired to realize quantum computing. A quantum gate 
is often realized by manipulating the parameters in the 
Hamiltonian of a system so that the time-evolution oper- 
ator results in a desired unitary gate. On the other hand, 
when the system has a degenerate energy eigenvalue, adi- 
abatic parameter control allows us to construct a quan- 
tum gate employing non-Abelian holonomy Jj] . Holon- 
omy corresponds to the difference between the initial and 
the final quantum states under an adiabatic change of 
parameters along a closed path (loop) in the parame- 
ter manifold M. Therefore, a desired quantum gate 
can be implemented by choosing a proper closed loop in 
Ai. This scheme is called the holonomic quantum com- 
puting (HQC). The idea was suggested first in Ref. Q 
and has been developed subsequently by many authors 
d, 0, 0, H, Hi ■ Holonomy is geometrical by nature, and 
hence it is independent of how fast the loop in the pa- 
rameter manifold is traversed. In addition, if the lowest 
eigenspace of the spectra is employed as a computational 
subspace, it is free from errors caused by spontaneous 
decay. Thus, HQC is expected to be robust against noise 
and decoherence [loj . 

In spite of its mathematical beauty, physical imple- 
mentation of HQC is far from trivial. The difficulties are 
(i) to find a quantum system in which the lowest energy 
eigenvalue is degenerate and (ii) to design a control which 
leaves the ground state degenerate as the loop is tra- 
versed. Several theoretical ideas have been proposed in 
linear optics 11 1, trapped ions [ll,[H,[l3|, and Josephson 
junction qubits [T5| . Recently, an experiment following 
the proposals made in [13, 13, where the coding space is 
not the lowest eigenspace, has been reported jig . 

Karimipour and Majd [l3| proposed HQC with a spin 



chain model. A single logical qubit is represented by two 
spin-i physical spins; the qubit is a dimer. In the dimer, 
the spins interact with each other through a Heisenberg- 
type interaction and its lowest eigenvalue is doubly de- 
generate for a particular set of the coupling constant and 
the Zeeman energies. The spin chain model may be ap- 
plicable to various physical systems, such as solid state 
systems. However, there is room to reconsider and im- 
prove their proposal. In particular, we want to inves- 
tigate whether the Heisenberg-type interaction is essen- 
tial for HQC in a spin chain model. Such considera- 
tion is necessary for clarifying the applicability of their 
proposal. In this paper, we closely follow the discussion 
given in Ref. [l7|| and construct holonomic quantum gates 
using isospectral deformations of an Ising model, instead 
of a Heisenberg model. In addition, we explicitly require 
that a path in the parameter manifold be closed for 
a holonomy to be well-defined. 

The paper is organized as follows. We briefly review 
holonomy associated with adiabatic time evolution of a 
quantum system in Sec. [Til We show how to construct 
one- and two-qubit gates as holonomies in Sections IIIII 
and IIVI Section |V] is devoted to summary. 



II. HOLONOMY 

Let us consider a Hamiltonian H acting on a Hilbert 
space H. (dimH — N), whose Ith eigenvalue is denoted 
as El. In particular, I = refers to the lowest eigenvalue. 
The Ith eigenvalue is gi-iold degenerate and its eigenvec- 
tors are written as \l, i) {i — 1, 2, . . . ,gi). We assume 
{I, i\m, j) = SimSij. Note that N 9i- 

An isospectral deformation of H is accomplished by 



H{t) = g{T)H g\T) (0 < r < 1), 



(1) 
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where (?(r) G U(7V). As a result, no level crossing takes 
place during the Hamiltonian deformation. We normalize 
T e [0, 1] so that iJ(0) = H{1) = H. This condition 
implies that the curve in be closed. The symbol r in 
Eq. (Ill) is the normalized dimensionless time: t ~ t/T 
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{0 < t < T), where T is the total time to traverse the 
loop. Note that T is long enou gh so that the adiabatic 
approximation may be justified |18| . 

We consider the particular isospectral deformation 

H{t) ^ e^^He-^\ (2) 

following [1, [l3| , where X is a constant anti-Hermitian 
matrix. We require [H, X] ^ and 

- a, (3) 

where 1 is the unit matrix. The first condition is required 
to implement non-trivial gates while Eq. ([3]) ensures that 
H{1) = H{0) = H. 

We readily obtain the instantaneous eigenvalues Ei (t) 
and eigenvectors |^ j;T) of H{t) as follows [sl. [iTj: 

Ei{t)=Eu |/, z; r) =e^n^ (4) 

which implies that no level crossing occurs during the 
deformation: Ei{t) ^ Ev{t) {I ^ r,Vr G [0,1]). We 
will exclusively work with the ground state {I = 0) from 
now on and drop the index I = whenever it cause no 
confusion. 

We use the unit in which h = 1 through the paper. Ac- 
cording to the adiabatic theorem, when the initial state 
\iP{t = 0)) is in the lowest eigenspace, the final state 
\tp{T = 1)) remains within this subspace. Then, we find 

|^(T = l)>=e-^^»^r|^(r = 0)), (5) 

where e"*^""^ is the dynamical phase and F e U(.go)- 
Equation ^ can be regarded as a gate operation con- 
necting the initial and final states. The unitary matrix 
r in Eq. ^ is the holonomy associated with the cyclic 
deformation Eq. ^ of the Hamiltonian and we write it 
as 

r = er^. (6) 

The anti-Hermitian connection A in Eq. ([6]) is given, in 
terms of X, as 

90 

A = ^(^|Xb•)|^)0■|. (7) 

In order to calculate A, the instantaneous eigenvector 
|0, i; t) in Eq. @ is substituted into 

A^Ar) = {^r\^\3■,r), (8) 

which is the definition of the ij component of the con- 
nection A. 

III. ONE-QUBIT GATES 
A. Hamiltonian 

We introduce a Hamiltonian 

i/lD = -UJCFiz - ^02z + ■h<^\zC!2z (9) 



as H in Eq. ([2]), where Oka is the a-component of Pauli 
matrices of the fcth spin (fc = 1, 2 and a = x,y, z). Equa- 
tion ([9]) corresponds to the Hamiltonian of two homoge- 
neous spin-^ particles interacting with each other via an 
Ising-type interaction as depicted in Fig. [TJ The strength 
of the interaction is parameterized by Ji, while that of 
the field by uj. We assume uj > and Ji > without loss 
of generality. 
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FIG. 1: Dimer consists of two spin-i particles. 

We denote the eigenvectors of dz as follows: crz\+) = 
\+) and crz\—) — —\—)- Furthermore, we introduce the 
following vectors which are eigenvectors of iJiD: = 
I + +>, iTo) = -i=(| + -) + I - +)), |r_) = I - -), and 

|5'o) = -ijd + -)-[-+)), where denotes |+)(»|+), 

and so on. The eigenvalues of |T+), |To), |T_), and |5o) 
are — 2ijj + Ji, — Ji, 2^ -I- Ji and — Ji, respectively. 

In case lj — Ji, there exists 3-fold degenerate lowest 
energy eigenvalue — Ji . The lowest eigenspace is spanned 
by the three eigenvectors |T+), |To) and \ So)- The unique 
excited state is |T_) and the energy difference between 
|T_) and the ground state is 4Ji. 

In this paper, we choose |0)l = 17"+) and ~ |Tb) 
as basis vectors of the logical qubit. Namely, the coding 
space for a single qubit is Ci = Span{|T+), |To)}. We 
note here that a qutrit may be implemented by using the 
three ground state eigenvectors. Although this coding 
is potentially interesting, it is beyond the scope of this 
paper. 



B. Implementation of one-qubit gates 

One-qubit gates are implemented by choosing X in 
Eq. @ as [13 

X = inn- {(Ti+(T2), (10) 

where n is a unit vector in M'^, while is a positive 
real number. It should be emphasized that undesired 
transitions into the irrelevant subspace Span{|S'o)} do 
not occur for this choice of X. This is easily seen from 
the identity (ctiq -I- cr2a)|'5'o) = 0. We obtain the anti- 
Hermitian connection restricted to the coding space Ci 
as follows: 

A\^^ = in[nz{lL + <JLz) + V2{nxcrhx + Uya-Ly)]. (11) 

where (Tl. = |r+)(To| + \To){T+\, aLy = -i{\T+){To\ - 
\To){T+\), = |T+>(r+|-|ro)(ro| and/L = |T+)(T+| + 
|ro)(2o|. Note that the presence of a closed loop in A4 is 
essential to obtain Eq. (jlip. 
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Using Eqs. ^ and (fTTj) . we find the unitary operator 
p — g-i^rij g-iO[nj(TLj+\/2(ri^c7L^+n„CTLy)] The Condition 
^^id(1) = HiY){0) is necessary for the closure of a loop 
in [19]. As for Eq. (fTO)) . this condition is obviously 
satisfied by taking f2 = ktt (k S N), where k is interpreted 
as the winding number of a loop in A4. 

Summarizing the above arguments, we find that X 
given in Eq. (jlOp implements a single qubit holonomic 
gate 

where m = , ^ , (^/2nx, \/2ny, Uz) and 6^ — 

a/2 — ^ y / 

K7r^2 — n^. In order that the holonomic quantum gate is 
non-trivial, the condition [i?iD, -'^l 7^ must be satisfied. 
This condition is equivalent to \nz\ 7^ 1. 

Note that the set of the rotation angles 0^ is discrete 
when Uz is fixed, and hence it is impossible to vary 0^ 
continuously in Eq. Nevertheless, we can find a 

quantum gate that approximates the de sired one with 
arbitrary accuracy provided that ^2 — n| is an irrational 
number [20| . 



C. Examples 



2. Arbitrary elements of SU(2) 

The holonomy T^i^'q Q-jC'*) is a quantum gate generat- 
ing a rotation around the x-axis by an angle 9^ — v^ktt. 
Although the set of 6^, is discrete, we can find k which sat- 
isfies \{9~9i^) mod 27r| < e for arbitrary 6 and e [2^1 . Fig- 
ure[3]shows the points (cos6'„, sin6'„) for k = 0, 1, . . . , 10. 
The point (1, 0) corresponds to k = (i.e., X = 0), in 
which no gate operation is performed. Similarly, we ob- 
tain an approximate rotation around the y-axis with ar- 
bitrary accuracy by taking m = (0, 1, 0). Consequently, 
it is possible to implement an arbitrary element of SU(2). 
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1. Hadamard gate 

The Hadamard gate, up to an irrelevant overall phase, 
is implemented by taking both | sin 6*^ | = 1 and m — 

(l,0,l)/\/2 in r^'(K). Note that Uz = ^2/3 and 6*^ = 
2K7r/\/3 for the choice of m and | sin6'K | — 1 cannot be 
satisfied exactly for any k. However we show that the 
Hadamard gate may be implemented with good accuracy 
for some k. Figure [2] shows sin^^ as a function of k, 
from which we find that | sin 0^1 — 1 for k = 3,10 and 
16. It can be proved easily that the Hadamard gate can 
be implemented with arbitrary accuracy by choosing a 
proper k. 
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FIG. 2: sin 6'k as a function of k. 



FIG. 3: Rotation angles 6k ~ 2ktv/\^ in F^J^'g o)(^) plot- 
ted as points (cos sin on the unit circle. The numbers 
in the figure indicate the values of n. 



IV. TWO-QUBIT GATES 

A. Hamiltonian 

A logical two-qubit system consists of two dimers. The 
Hamiltonian is 

H2D = H' + H^, (13) 

where = —Jictiz — Ji cr2z + Ji<yiz^2z and = 
— J20'3z — ^2 CT4z -|- J20'3zO'4:z- Equation ([15)) is used as 
H in Eq. (g]). Here {H'^) is the Hamiltonian of a 
single dimer to which the first and the second (the third 
and the fourth) spins belong. We easily find the 9-fold 
degenerate lowest eigenvalue — Ji — J2. 

We take the following coding space C2 
for the logical two-qubit system: C2 = 
Span{|T+)i|T+)2, |r+)i|To>2, |To)i|r+)2, |To)i|To)2}. 
The vector |r+)i denotes the eigenvector | + +) 
associated with H^, for example. 
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B. Controlled-e'"^ gate 

We choose the fonowing generator of the isospectral 
deformation in Eq. ([2]); 



(14) 



where ~ iriiVli ■ (cri + cr2), X^ = in2^2 ■ (C3+C4) and 
X^-'^ = iJ{ai^a3z + (Ji^aAz + 027.0^^ + cf2zOaz)- Here, 
ni and n2 are unit vectors in M^, while r^i, 172, and J 
are positive real numbers. No undesired transitions into 
the non-coding space take place for this choice of X. We 
elaborate this point in Appendix IB] 

Let us find the corresponding anti-Hermitian connec- 
tion with the following unit vectors n\ and n2; 



m = (0, 0, 1 
n2 = 



1 - ni, 



0, ^22 



"22 



J 



(15) 
(16) 

(17) 



We assume Q.2 > J, which guarantees [H2D, X] 7^ 0. 
Then, we obtain the following anti-Hermitian connection: 



V2^2n2xlL ® CLx + Jc^hz ® <y-Lz] ■ (18) 



Note here that ^2n2x — \/^\ — J^- 

The condition e-^ — restricts the parameters as 



O2 = n+TT (k+ e N), 



(19) 
(20) 
(21) 



8J2 ^ K_7r (k_ G N), 
Vli = k'tt [k' e N). 

Their derivation is given in Appendix [B) It follows 
from Eqs. dH]) and ([20]) that k_ > k+ and J = 

Y^K?_ — k;+. Note that the assumption 1^2 > is 

equivalent to 3k+ > k_. 

As a result, we obtain a two-qubit gate 

r(2)(K:, k') = (-i)«'rLU(K;, K')rC(K), (22) 

where k = {k+, k_). The local unitary gate T^^ {k, k') 
is defined as 



r 

where v 



LU 



(/«, k') = e-*(«''^+^)'^i-- 



(23) 



= A/2K^7r2 — J2 and vk — 
(\/2K^7r2 — 2 J2, 0, J). The essential nonlocal oper- 
ation is r^(/«), which is a controlled-e'^^ gate and is 
given by 



|i)t,(i| 



i2J<TLj 



(24) 



Thus, r(2)(/t, k') followed by the local unitary gate 
T^^{k, n'y implements the controUed-e*^'^ gate with 
9 = 2J 
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FIG. 4: Rotation angle 2 J in the controlled-e'"^, T^Ik). (a) 
The value of 2J with the unit of tt is plotted with respect to 
K_, when the fixed values of k+ are given, (b) The values of 
cos 2 J and sin 2 J are plotted as the points on the unit circle. 
The plotting marks are the same meaning as in Fig.|4l[a). 



The rotation angle 2 J in r^(K) is characterized by 
and K_. Wc can, however, find and k_ which 
satisfies 1(6* — 2 J) mod 27r| < e for arbitrary 9 and e [23 |. 
We show the attainable values of 2 J for the various sets 
of (, 

be chosen such that k+ < 



K_) in Fig. m (a). Note that (k+, k„) have to 
< 3K4. is satisfied. We 
calculate (cos 2J, sin 2 J) as in Fig.|3](b). Thus, a tunable 
coupling control scheme, in which J is controllable, is 
necessary to achieve various rotation angles in r^(/«). 

Alternatively, by repeating the above gate n times, we 
can implement the controlled-e*^"'^'^^ gate. For an ir- 
rational J and a given 6, it is possible to find n such 
that \{2nJ — 9) mod27r| < e for an arbitrary e. In this 
way, one may implement the controUed-e**'^^ gate with 
arbitrary precision. In particular, we can implement the 
controUed-Z gate, which is a constituent of the universal 
set of quantum gates, by taking 9 — Tr/2. 



V. SUMMARY 

We have constructed holonomic quantum gates using 
isospectral deformations of an Ising model Hamiltonian. 
These gates are the Hadamard gate, rotations around x- 
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and y-axes and the controlled-e'^'^ gate. The closure of 
the loop in the parameter manifold leads to a discrete 
set of gates. These gates are, however, dense in SU(2) 
and SU(4) and any one- and two-qubit gates can be im- 
plemented with arbitrary accuracy. A spin-chain model 
is a good candidate to implement holonomic quantum 
gates with the ground state eigenspace. We will propose 
a more feasible scheme based on the present analyses in 
our future work. 



Accordingly, the time evolution operator at t = 1 is 

;7(T-l) = e^e-*(-^^+^^). (A3) 



When the adiabatic approximation is valid and the initial 
state is in the lowest eigenspace, we obtain 
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APPENDIX A: EMULATION OF HOLONOMIC 
QUANTUM GATES 

An exact solution of the Schrodinger equation is ob- 
tained when the Hamiltonian is deformed according to 
Eq. ([2]). Let us consider the Schrodinger equation in term 
of the dimensionless time r: 

^-^mr))^THirmr)), (Al) 

where H{t) = e^'^ H e~-^'^ . Introducing \4>{t)) = 
e~'^'^\ip{T)) , we obtain the following equation for |0(t)): 

z^mr))^{-iX + H)\<Pir)). (A2) 



f/(T = l)Po - e-'^«^rFo, (A4) 



where Pq is the projection operator on the lowest 
eigenspace. From Eqs. (|A3|) and (|A4p . we obtain 



e''e-'^-'''+''^^Po = e-'^"^TPo. (A5) 



The right hand side of Eq. (|A5p amounts to a holonomic 
quantum gates, up to the dynamical phase e^'^°-^. The 
left hand side of Eq. (|A5|) may be interpreted as imple- 
mentation of the quantum dynamics by a pulse sequence. 
Equation. (jASp provides a way how to emulate a holo- 
nomic gate r with a series of pulses of g-i-i-^^+^T) g^j-^^ 
e-^ , although it may not be a genuine realization. 
APPENDIX B: GENERATOR OF ISOSPECTRAL 
DEFORMATION FOR TWO-QUBIT GATES 

Let us analyze the the generator of the isospcctral de- 
formation for two-qubit gates in Sec. IIVBI In this Ap- 
pendix, we denote the eigenvector |T'+)i|r+)2 simply as 
|Tl|-T+), for example. 

First, we focus on and A^ in Eq. (fT4|) . We obtain 



2niz V^{nix - iniy) 

= iVti I ^{m^+inxy) 





2n2z \f2{n2x - in2y) 

X'^ = I V2in2^ + in2y) 





(Bl) 



The matrix element {X^)ki corresponds to i(fc|A'|^)i, 
where |1), = \T+),, |2), = |To)„ and |3), = |5o). 
(i = 1, 2). One can easily find that {T+So\X^\T+Tq) = 0, 
for example. This means that undesired transitions into 
the subspace irrelevant to quantum computation never 
occur under X^ and A^. 

Next, let us consider A^"^ in Eq. p4| . The following 
observations are useful to calculate the matrix elements: 

ak.\T+), = |T+)„ (B2) 
<TkATo), = {~\f+^\Sa),, (B3) 



where /c = 1, 2 for i = 1 and fc = 3, 4 for i — 2. 

Thus, we only have to consider five ele- 
ments {T+T+\X^-^\T+T+) , {T+S^\X^--'\T+Tq) , 

(5or+|Ai-2|ror+), {S^S^\X^-^\nSo) , and 

(5oro|A^"^|To5o). Note that one can easily calculate 
the other non-trivial elements (e.g., {ToT+\X'^-'^\SoT+)) 
from the anti-Hermitian property of A^"^. The other 
elements which are related to the lowest eigenspace 
of H2D trivially vanish. As a result, we obtain 
the following results: {T+T+\X^-'^\T+T+) = i4J 
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and (T+^ol^'-'lT+To) - {SoT+\X^-^\ToT+) = 
{SoSo\X^-^\ToSo) = {SoTo\X^-^\ToSo) =0. As a 
result, all the matrix elements corresponding to the 
undesired transition to the irrelevant subspace vanish. 

Finally, we calculate when the unit vectors rii and 
n2 are given by Eqs. pS]) and (fT6|) . respectively. We ob- 
tain 

- |r+)n(r+|®e^^+ + |T_)n(r-|®e'^- 

+ (|ro>ii(ro| + |5o>ii(5o|)®e^^°, (B4) 
Yo ^ f^ana • (0-3 + 0-4), (B5) 
Y± = ±217il2 +iy±k±- (0-3 + 0-4), (B6) 



where i/+ = ^12,1^- = vf^l + SJ^, and I2 = 
|r+)22(T+| + |To)22(ro| + |T_)22(r_| + \So)22{So\- In 
particular, the values of i'± is important for imposing 
the condition = 1. The unit vector k± is given 
by h'±k± = {^2n2x, 0, 02n2z ± 2J). The condition 
= 1 impHes e*'*'" = 12 and e''*'± = 12. The con- 
dition e*^" = IL2 implies there is k+ G N such that 
CI2 — k+tt; this requirement is equivalent to Eq. (|19p . 
Similarly, if Eqs. ((20l) and pT|) are satisfied, then we ob- 
tain e^'*'± = la. 
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